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Abstract: The measurement of Baryonic Acoustic Oscillations from galaxy surveys is 
well known to be a robust and powerful tool to constrain dark energy. This method 
relies on the knowledge of the size of the acoustic horizon at radiation drag derived from 
Cosmic Microwave Background Anisotropy measurements. In this paper we quantify the 
effect of non-standard initial conditions in the form of an isocurvature component on the 
determination of dark energy parameters from future BAO surveys. In particular, if there 
is an isocurvature component (at a level still allowed by present data) but it is ignored 
in the CMB analysis, the sound horizon and cosmological parameters determination is 
biased, and, as a consequence, future surveys may incorrectly suggest deviations from a 
cosmological constant. In order to recover an unbiased determination of the sound horizon 
and dark energy parameters, a component of isocurvature perturbations must be included 
in the model when analyzing CMB data. Fortunately, doing so does not increase parameter 
errors significantly. 
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1. Introduction 

Born on the basis of simplicity, the adiabatic framework describing the initial conditions 
for the perturbations, consolidated its popularity when it was shown to be predicted by 
the simplest one-field inflationary model [1, 2]. Although this simplest and basic adiabatic 
picture is widely accepted and provides an excellent fit to current data (e.g., [3]), there is 
not a priori reason to discard different and more general initial conditions involving entropy 
isocurvature perturbations. Such models also relay on motivated physical assumptions, e.g. 
models of inflation with more than one field [4, 5, 6, 8, 9, 7, 10], neutrinos isocurvature 
perturbations [11], axion dark matter [12, 13, 14, 15, 16, 17, 18, 19, 20] or the curvaton 
scenario [21, 22, 23]. 

Even if pure isocurvature models have been ruled out [24, 27, 25, 26], current obser- 
vations allow for mixed adiabatic and isocurvature contributions [3, 28, 29, 30, 31, 32]. As 
shown in [36, 30, 34, 37, 35, 38, 39], the initial conditions issue is a very delicate prob- 
lem: in fact, for current cosmological data, relaxing the assumption of adiabaticity reduces 
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our ability to do precision cosmology since it compromises the accuracy of parameter con- 
straints. Generally, allowing for isocurvature modes introduces new degenerations in the 
parameters space which weaken considerably constraints. 

The Cosmic Microwave Background radiation (CMB), being our window on the early 
universe, is the preferred data set to learn about initial conditions. Up to now, however, 
the CMB temperature power spectrum alone, which is the CMB observable better con- 
strained so far, has not been able to break the degeneracy between the nature of initial 
perturbations (i.e. the amount and properties of an isocurvature component) and cosmo- 
logical parameters, e.g. [34, 35, 40, 41]. Even if the precision measurement of the CMB 
first acoustic peak at ^ ~ 220 [25, 26] ruled out the possibility of a dominant isocurvature 
mode, allowing for isocurvature perturbations together with the adiabtic ones introduce 
additional degeneracies in the interpretation of the CMB data that current experiments 
could not break. Adding external data sets somewhat alleviates the issue for some degen- 
eracy directions e.g. [28, 31]. As shown in [42], the precision polarization measurement 
of the next CMB experiments like Planck will be crucial to lift such degeneracies, i.e., to 
distinguish the effect of the isocurvature modes from those due to the variations of the 
cosmological parameters. 

It is important to keep in mind that analyzing the CMB data with the prior assumption 
of purely adiabatic initial conditions when the real universe contains even a small isocurva- 
ture contribution, could lead to an incorrect determination of the cosmological parameters. 
In fact, the presence of an isocurvature component changes the shape and the location of 
the CMB acoustic peaks, mimicking the effect of parameters such as ^mh^, Hq and w (see 
[30] and references therein). 

Here we investigate wether such an effect has an impact on "standard rulers" , like the 
sound horizon at recombination, inferred from CMB observations. 

This issue turns out to be very delicate and important especially in view of the next 
generation of galaxy surveys which aims at probing with high accuracy the late time expan- 
sion and thus the nature of dark energy by means of Baryon Acoustic Oscillations (BAO) 
at low redshift (z < 2). 

The BAO in the primordial photon-baryon fluid leave in fact an imprint in the large 
scale matter distribution and, at each redshift, their physical properties are related to the 
size of the sound horizon at radiation drag (i.e., when baryons were released from the 
photons) e.g. [43] . Thus it is possible to use measurements of these oscillations at different 
redshifts as standard rulers. The key idea is that measuring the BAO signal at low redshift 
with large scale structure surveys and knowing the size of the sound horizon at radiation 
drag with high precision from the CMB should allow to probe the expansion history, i.e., 
the Hubble parameter H{z), and the angular diameter distance Da{z) at different redshift 
[44] and thus the dark energy properties. 

It has been found that the effects of a number of theoretical systematics such as non- 
linearities, bias etc, on the determination of the BAO location can be minimized to the 
point that BAO are one of the key obscrvablcs of the next generation dark energy ex- 
periments. It is therefore important to investigate the robustness of the method to other 
theoretical uncertainties. A crucial assumption is the possibility to measure accurately and 
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robustly the acoustic sound horizon at radiation drag. Rcf [45] showed that generally CMB 
observations yield a robust determination of such standard ruler although there could be 
possible systematic effects introduced by deviations from the standard evolution of the 
early universe. Refs [46, 47] have considered two early universe non-standard phenomena 
that change the standard recombination process and could affect BAO interpretation. In 
this paper we investigate another possible deviation from the minimal, standard cosmo- 
logical model that could affect the BAO interpretation as standard rulers if not taken into 
account: the presence of an isocurvature contribution to the initial conditions for primor- 
dial perturbations. This point constitutes the new part of this work, since such an effect 
has not been explored or quantified in details so far. 

The paper is organized as follows: in the next section we give a brief introduction to 
isocurvature modes and the notation used. This is only review material but it is reported 
here to clarify our set-up and notation. In section 3 we investigate the effects on forecasted 
parameter constraints for an experiment with characteristic similar to the Planck satellite, 
arising from the presence of an isocurvature contribution to the primordial perturbations. 
(In what follows we will refer to such an experiment as "Planck" for short). We quantify the 
systematic effects on quantities like the sound horizon at radiation drag and the expansion 
history parameter, introduced if an isocurvature component at a level allowed by current 
data, was present in the data but ignored in the analysis. 

In section 4 we present our results on the impact of isocurvature contributions on the 
BAO observables and we explore the implications on dark energy constraints from future 
surveys forecasts. Finally in section 5 we present our conclusions. Appendix A reports the 
relevant equations and definitions for BAO, in Appendix B we recap the basic equations 
for the Fisher matrix approach to forecasting and Appendix C explores the source of CMB 
parameters degeneracies introduced by an extra parameter for the isocurvature contribution 
to the initial perturbations. 

2. Isocurvature: some theory and notation 

The simplest characterization of the primordial perturbations ^ is described by the adiabatic 
framework: all particles species are perturbed in spatially uniform ratio so that the relative 
ratios in the number density remain unperturbed. Given two different particle species, X 
and Y, the adiabatic condition requires: 



where nx and ny are the correspondent particle number densities. 

The global perturbation on the matter component, because of the Einstein equations, 
is associated with a curvature perturbation. This is the reason why adiabatic perturbations 
are also called curvature perturbations. 

^In this context primordial perturbations refers to perturbations defined deep in the radiation era on 

superhorizon wavelengths with modes initially excited well before recombination so that any decay mode 
had time to vanish and the main cosmological components are: photons, baryons, neutrinos and CDM. 




(2.1) 
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However, this is not the only possibiHty. There can exist in fact perturbations associ- 
ated with fluctuations in number density between different components of the cosmological 
plasma in the early universe, well before photon decoupling. These are called isocurvature 
or entropy perturbations and are generated by the stress fluctuations through the causal 
redistribution of matter under energy-momentum conservation. Causality and momen- 
tum conservation require that the initial curvature perturbations must vanish (spatially 
varying abundance of particles species are arranged to cancel locally leaving the curvature 
of the spatial hypersurface unperturbed). Density perturbations are then produced by 
non-adiabatic pressure (entropy) perturbations. See [48] for a detailed description. 

The entropy perturbation between two particle species (i.e., fluid components) can be 
written in terms of the density contrast 5i = ^ and the equation of state parameter Wi as 

Pi 

2. 

SxY = - (2.2) 

1 + WX 1 + WY 

Formally, this is equivalent to: Sxy = — which quantifies the variation in the 
particle number densities between two different species. 

More precisely (see [11] and references therein) the most general description of a pri- 
mordial perturbation accounts for 5 non-decaying (regular) modes corresponding to each 
wavenumber: an adiabatic (AD) growing mode, a baryon isocurvature mode (BI), a cold 
dark matter isocurvature mode (GDI), a neutrino density (NID) and a neutrino velocity 
(NIV) mode. For the dark matter component ^ = so that the CDM isocurvature 
mode can be written as: ^ 

Sc = dc- -^7, (2.3) 

where 6x = is the energy density contrast of the X particle species. Analogously the 
baryon (b) and the neutrinos (i^) isocurvature modes take the form: 

S, = S,- ^S^ (2.4) 

5. = - Is^. (2.5) 

Of course for the adiabatic mode Sc = S,^ = Sh = 0. 

One of the crucial distinctions between the adiabatic and the isocurvature models relies 
on the behavior of the fluctuations at very early time during horizon crossing at the epoch 
of inflation (or analogous model for the very early universe) . In the adiabatic case constant 
density perturbations are present initially and imply a constant curvature on super-horizon 
scales, while in the pure isocurvature case there are not initial density fluctuations which 
instead are created from stresses in the radiation-matter component. 

2.1 A worked example: the curvaton model 

An alternative to the standard single field inflationary model is the curvaton scenario 
[21, 22, 23]. This model relies on the inflaton, the light scalar field that dominates the 



^Recall that the continuity equation, which follows from the energy-momentum conservation, takes the 
form p = -3H{P + p) = -3Hp{l + w). 
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background density during inflation, and the curvaton field a which, decaying after infia- 
tion, seeds the observable cosmological perturbations. Within this framework, the initial 
conditions then correspond to purely entropy primordial fluctuations, namely isocurvature 
perturbations, because the curvaton field practically does not contribute to perturbations 
of the metric. If the field a decays sufficiently early and transforms completely into ther- 
malized radiation, pure adiabatic perturbations with a blue spectrum are generated and 
the initial isocurvature density perturbations disappear. 

In general, by tuning the decay dynamics of a, the curvaton scenario allows for mixed 
adiabatic and isocurvature fluctuations, with any residual isocurvature perturbation cor- 
related or anti-correlated to the adiabatic density one and with the same tilt for the both 
spectra: Had = niso- 

Besides the fact that current data are compatible with this theoretical picture, the 
curvaton model attracted growing attention because it predicts primordial non-gaussianity 
features of the local type in the spectrum of primordial perturbations. In the paper we 
use the curvaton scenario as a working example for a model that gives rise to a small 
fraction of correlated GDI isocurvature. As shown in the next sections, our analysis can be 
generalized to models with an arbitrary amount and type of isocurvature. The curvaton 
has Uad = niso, therefore our findings are quantitative only for this case: we will comment 
on how our findings should qualitatively hold also for Uad 7^ niso- 

2.2 Isocurvature and the cosmological parameter fiso 

A common parametrization for the isocurvature perturbations is given by [49] : 

defined as the ratio between the entropy S and the curvature (adiabatic) TZ perturbations 
evaluated during the radiation epoch and at a pivot scale ko. In our case we set 0.05 Mpc~^, 
as often done in the literature. The correlation coefficient can be then defined in terms of 
an angle Ajt^ such that: 

Throughout this paper we will use a sign convention such that fiso < will correspond 
to correlated modes. As pointed out in [34], in a practical implementation of a multi- 
parametric analysis using Markov Chains Monte Carlo (MCMC), the fiso parametrization 
(instead of favors small multiplier in front of the isourvature component thus adding 
an implicit bias towards pure adiabatic model in the posterior. Another parametrization in 
terms of the amplitude parameter a for which the data have a linear response is given by: 
a = (^i^j^i y P = cos(Afcp), where maximally correlated (anticorrelated) modes correspond 
to P = +1 { P = —1). Our analysis is mostly concerned with degeneracies directions and 
is thus less sensitive to the effect of the prior choice on the posteriors. 

For the curvaton scenario with mixed adiabatic and purely correlated CDM isocurva- 
ture modes which we chose as working example, the most recent constraints are a-i < 0.011 
(95% CL), [3]. 
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Figure 1: The different components of tlie unit- amplitude temperature angular power spectrum 
defined in Eq. 2.11 (dot-daslied line) for the cosmological parameters listed in table 1 and in the 
case TiigQ — Tifid : Cf^ (solid black line), C/^^ (dashed hue) and C^°'' ( orange line). See Eqs. (2.8 
- 2.10). 



Since we are interested in the effect of adding an isocurvature contribution on the CMB 
angular power spectrum it is useful to give the explicit expression for this. In general the 
two-point correlation function or power spectra for the adiabatic mode, the isocurvature 
mode and their cross-correlation can be described by two amplitudes, one correlation angle 
and three independent spectral indexes {uad, niso and ricor), so that, in the case of the 
CMB, the respective power spectra are: 



cr=/y[0mf(^) , (2. 



dk r^io,, / , / k 



and 



icor I dk ^ ^ I k 



C^or ^ J -ea.(^)0.so(^) ^__J _ (2.10) 

Here &'^'^{k) and 0£*°(A:) are the radiation transfer functions for adiabatic and isocurva- 
ture perturbations that describe how an initial perturbation evolved to a temperature or 
polarization anisotropy multipole i. The total angular power spectrum takes then the form: 

Ci = {nl,)[Cf + fLcr + 2^oCos Afc^Cn. (2.11) 

Recall that in this paper we will always work in the case where Uad = niso- In Fig- 1 we 
show the unit-amplitde components of the temperature angular power spectrum listed in 
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the above equations: Cf (solid black line), Cj (dashed line) and C^°'^ (orange line). 
The cosmological parameters used are given in tab. 1 and in the total mixed ad+iso power 
spectrum (dot-dashed line) the isocurvature fraction is fiso = —0.08. Prom the plot and 
from eq. (2.11) it is clear that, for small isocurvature fractions, the main isocurvature 
contribution come from the mixing term coefficient 2/jsoCOS A^g. 

3. Analysis 

The aim of this section is to investigate and quantify the following. Assume that the Uni- 
verse contained a small isocurvature contribution to the primordial perturbation at a level 

allowed by current constraints. If forthcoming Planck satellite data were to be analyzed in 
the context of purely adiabatic initial conditions: a) what would be the systematic error 
on quantities like the sound horizon at radiation drag and what would be the implications 
for the reconstructed expansion history from BAO observations? (sec. 3.1) b) what would 
be the shift on cosmological parameters and on the sound horizon at radiation drag as a 
function of the amount of isocurvature contribution? How does this shift compare with the 
forecasted statistical errors on the same quantities? (sec. 3.2). 

Throughout we will assume a flat universe. Recall however, as pointed out in [30, 31], 
that assuming spatial flatness in isocurvature studies strongly biases the result toward pure 
adiabaticity. 

3.1 MCMC approach to forecasts for a Planck- like experiment 

We adopt a modified version of the CAMB code [33, 28, 29] which includes the initial 
conditions for the correlated mixed adiabatic (AD) and cold dark matter isocurvature 
(GDI) modes predicted by the curvaton scenario. With this code we generate two fiducial 
models: an adiabatic one and one with an isocurvature contribution still allowed by current 
data, as summarized in Table 1: 

• Fid.AD- Fiducial model Adiabatic (AD) - 7 parameters. 

• Fid.ISO- Fiducial model curvaton - 8 parameters. 

The set of parameters we consider is given by: = flb^^, uic = ^ch"^ that are respectively 
the baryon and cold dark matter physical density fractions, the optical depth at reionization 
T, the Hubble parameter Hq = 100 hKm sec Mpc~^, the scalar amplitude Ag = ln[10^'^Tlrad] 
, being TZrad the curvature perturbation during the radiation era, the scalar adiabatic tilt 
Us, the dark energy equation of state parameter w and fiso the isocurvature fraction. 

We estimate Planck satellite errors on the CMB temperature and polarization and the 
form for the likelihood using the specifications of Ref. [51, 52] and the technique presented 
in [53]. 

We use the code cosnioMC [56] to sample the posterior distribution for the cosmological 
parameters and a number of derived parameters using the Markov Chain Monte Carlo 
(MCMC) method. 

Each fiducial is sampled with two different models: 
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Table 1: Adiabatic and isocurvature fiducial models parameters. 



Parameter 


Ad-ACDM 


Isocurvature 




U.UZZ 


U.UzZ 


OJc 


0.11 


0.11 


h 


0.704 


0.704 


T 


0.073 


0.073 


Us 


0.96 


0.96 


As 


2.410"^ 


2.410"^ 


w 


-1 


-1 


fiso 





-0.08 



• fit.AD Model with parameters: {uJb,cOc,0,TjiT'STAs, {w)} {fiso = 0) 

• fit.ISO Model with params: {coh, coc, t, Us, As, (w), fiso}- 

A flat prior on these parameters is assumed. As usual the variable 6 = lOO^^^^y is used 
(and a fiat prior on this variable is assumed) instead of Hq or fl/^ since it is a "physical" 
parameter and improves MCMC performance; rs{z*) and Da{z*) are, respectively, the 
sound horizon and the angular diameter distance at decoupling (see Eqs A.3 and A.l). 
Both the AD and ISO analyses are run once fixing w = —1 and once letting w ^ —I (but 
still keeping it constant in time). 

For our purpose we will mainly use the chains: i) fiducial model curvaton fitted with 
adiabatic (Fid.ISO-fit.AD) for both cases cosmological constant {w = —1) and w ^ —1 but 
constant (Fid.ISO-fit.ADiu) and b) fiducial model adiabatic fitted with an anti-correlated 
isocurvature model (Fid.AD-fit.ISO) in the w = —1 case. (The other cases Fid. ISO-fit. ISO 
and Fid.AD-fit.AD are used only for cross checks and for testing that our procedure is not 
biased). The recovered best fit values (N-dimensional maximum likelihood location) for 
the cosmological parameters are listed in table 3.1. Only in the case where the mean best 
fit model is different from the maximum likelihood one both are reported. 

The "goodness of fit" parameter is 0.193871 and 0.464991, respectively for Fid.ISO- 
fit.AD and Fid.AD-fit.ISO assuming that Planck has 4000 independent data points (as- 
sumed to be the for temperature and polarization) so that both fits are acceptable. 

To begin with, we are interested in disentangling the effects of an (ignored) isocurva- 
ture contribution on an 'early' type observable like the sound horizon at the radiation drag 
and on a 'late' type observable like the expansion history at low redshift {0 < z <2). The 
acoustic oscillations imprinted in the galaxy distribution can in principle be measured at 
different redshifts so it is important to distinguish the effects of an isocurvature contribu- 
tion on expansion-history observables (e.g., H(z)) or determining the initial layout of the 
acoustic feature {rs{Z(i))- 

^Q{i^,X^) = 1 — r(j'/2,x^), where u — n — m are the degrees of freedom (the difference between the 
number of independent data points n and the number of parameters m) and is the minimum x^- 
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Fid.ISO-nt.AD 


Fid.ISO-nt.ADt/; 

best lit 


Fid.ISO-fit.ADiu 
mean 


Fid.AD-nt.ISO 


AlnL 


68.997 


67.784 




0.198 




0.02189 


0.0219 


0.02189 


0.022 


iOc 


0.iUo49 


O.IOODOO 


O.lOool 


nil nnoT 
0.110087 


Us 


0.9784 


0.97804 


0.9783 


0.9597 


w 


-1 


-1.342 


-1.171 


-1 




0.249 


0.171472 


0.2113 


0.26691 




0.757084 


0.828528 


0.7886 


0.733089 




72.4005 


86.2415 


79.52 


70.3461 



Table 2: MCMC-recovered maximum likelihood parameter values (i.e. likestats file). For the 
varying w chains we report the mean (.margestats file) parameters since in this case they differ 
from the maximum likelihood ones (in the other cases there is no noticeable difference), a: In this 
case the best fit model has a tiny anti-correlated isocurvature model; the best fit gives a slightly 
negative fiso but well within 1 — a error. 



3.1.1 Sound horizon at radiation drag rs{zd) 

The sound horizon at radiation drag rs{zd). i.e., the distance sound can travel up to the 
time when the baryons were released from the photons, is defined by Eq. (A. 3). The 
radiation drag redshift Zd, according to the fitting formula from [43] used in the literature 
so far, can be explicitly written as: 

1291(^)^fe^)°-^5i 
= l + 0.659(n^hT^^3 [l + bii^^h ) ^] (3.1) 
h = 0.313(n^/i2)-o.4i9[j ^ o.607(n^/i2)°-^^^] 

b2 = 0.238(O^/l2)0.223^ 

with Ujn = + ^cdm the matter density parameter. 

Actually this is an approximation and the drag epoch za can be calculated more pre- 
cisely, as pointed out in [58], as the epoch at which the drag optical depth equals one: 

Td{zd) = T— / dz- — = 1, (3.2) 

4 Jo da 1 + z 

where Xe{z) is the fraction of free electrons and ax the Thomson cross-section. 

By using the best fit values for the parameters for each type of chain (table 3.1), we 
calculate the sound horizon at radiation drag numerically in both cases. For the second 
more precise calculation we used a modified version of CAMB to extract both Zd and rs{zd) 
directly from the code. We find that, the use of the fitting formula yield values for rs{zd) 
inaccurate by up to 4%. The results are summarized in the table 3.1.1. It is clear from 
the 7th column of table 3.1.1 that if the underlying fiducial model is the curvaton and it 
is fitted with an adiabatic model the systematic shift in the estimated value of rs{zd) is 
~ 1% for both the w = —\ and w cases, while if the adiabatic fiducial model is fitted with 
an (anticorrelated) isocurvature model -as expected- there is no systematic shift. 
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Type of chain 




^CAMB 


^s[Z^ ) 




Fid. Adiabatic 


1018.7948 


1058.5975 


154.198 


150.449 


Fid.AD-fit.ISO 


1018.9217 


1058.559 


153.088 


150.430 


Fid. Curvaton 


1018.8711 


1058.52127 


154.198 


150.449 


Fid.ISO-fit.AD 


1018.1467 


1057.910919 


155.731 


151.888 


Fid.ISO-fit.ADu; best fit 


1018.7948 


1057.94906 


153.803 


151.823 



Table 3: Radiation drag redshift Zd and the sound horizon rs{zd)' we calculate the radiation drag 

redshift Zd by using the fitting formula from [43] and by using a modification of the CAMB code 
[58]: z^^'^ and z'^^'^^ . In the last two column we report the sound horizon at radiation drag rs{zd) 
(Mpc) for both cases. Note that (see sec. 3.2) the 1 — a error forecast for Planck on rs{zd) is 0.359 
(Mpc). 



Adding an isocurvature contribution implies a small variation in the estimation of the 
sound horizon at radiation drag that however should be taken into account for the Planck 
experiment since it is expected to measure such observable with high accuracy (better 
that 1%) , as we will discuss in §3.2. The implication of this systematic shift on the 
interpretation for future BAO data is explored in §4. 

3.1.2 Expansion history H{z) 

By using the best fit values for the parameters of each chain, we calculate the expansion 
history H{z) up to redshift 2; = 2, which is the redshift region of interest for studying 
the late-time expansion and accessible by current and future large-scale structure BAO 
surveys. Within this range and in the case of a flat model it takes the form: 

H{z) = Ho^nmil + zf + J7a(1 + z)3a+«'). (3.3) 

We find that the variation in the CMB-predictcd expansion history H{z) due to the 
fact that the underlying fiducial curvaton model is fitted with a pure adiabatic model is 
about 3% for a ACDM model at z = while it can reach 25 % for a model with varying 
dark energy equation of state w. 

These results are a quantitative confirmation of what previously anticipated: working 
with the prior assumption of purely adiabatic initial conditions when the real universe 
contains an isocurvature contribution could lead to an incorrect determination of the cos- 
mological parameters. In particular, the presence of an isocurvature fraction can affect the 
determination of the value of the sound horizon at the radiation drag. It also affects the 
CMB-recovered value of the Hubble parameter at the present time Hq and the expansion 
history -H{z)- inferred from CMB analysis. 

3.2 Planck Fisher analysis 

The MCMC approach is well suited to explore a single case (for a particular choice of 
fiso value). To study the dependence of the effect illustrated in §3.1 on the size of the 
isocurvature contribution, a Fisher matrix approach is better suited. The Fisher matrix 
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Channels (in GHz) 


44 


70 


100 


143 


217 


Beam FWHM 


24 


14 


10 


7.1 


5 


Final noise per arcminute fiK 


180.36 


180.74 


77.07 


57.46 


94.42 



Table 4: Planck experiment setting (upper table) and experimental specifications for the surveys 
used in this paper (lower table): redshift range z, survey area (A) in squared degrees, fraction of 
the sky fsky and mean galaxy number density n 





BOSS 


EUCLID 


z 


< 0.7 


< 2 


A {deg^) 


8000 


20000 


fsky 


0.2 


0.8 


n 


2.66 10-^ 


1.53 IQ-s 



analysis [54] (see B for more details) set up for the Planck experiment is calculated for a 
fiducial model with fiducial parameters: 

{r, ns,dns/d\n k, Zre, cob, oJc, h, Asjiso} = {0.01, 0.963, 0., 0.84, 0.02273, 0.1099, 0.72, 0.8169, -0.01}, 

(3.4) 

where Ug, dus/dln k and Zre are, respectively, the spectral index, the running of the spectral 
index and the rcionization redshift. 

In appendix B we review the the Fisher matrix error estimation procedure in the case 
of a CMB experiment. In the application here we use the modified CAMB version used in 
§3.1 to compute the Q. 

Beyond error estimations and forecasting, there exists another powerful tool encoded 
in the Fisher formalism. This enables one, without recomputing the full covariance matrix, 
to calculate the shift in the best fit parameters SOa due to the fact that an arbitrary number 
of parameters has been fixed to a wrong fiducial value. See [55] for more details. 

In general, given a number p of parameters ^'-y (7 = 1, ..,p) fixed to an incorrect value, 
which differs from the "true" value for an amount S^j, the resulting shift on the best fit 
value of the other n parameters 9a {a = 1, ...,n) is: 

Here is the sub-matrix of the inverse Fisher matrix, corresponding to the 6a 

parameters (i.e. the inverse of the Fisher matrix, without the rows and columns corre- 
sponding to the "incorrect" parameters) and S^^ is the Fisher sub-matrix including also 
the ^'-y parameters. 

Here in particular we are interested in computing the shift in the best fit parameters 
induced by setting the amount of isocurvature fiso to an incorrect value: S'^j = Sfiso and 
0a = {r,ns,dns/dlnk,Zre,uJb,u;c,h,As}. 

Moreover we want an estimate of this effect on the value of the sound horizon at 
radiation drag rs{zd), which is a parameter not directly included in the Fisher matrix. 
This can be done via a change of base in the Fisher approach. We calculate rs{zd) (see 
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Eq. (A. 3) numerically with values for the cosmological parameters given by Eq. (3.4). In 
general, if p = {pi, ...,pn) are the parameters in the Fisher matrix and q = (gi, qq) 
are the new ones, the new Fisher matrix F?^^^ respect to these will be: 



In our case the Jacobian matrix J„j = will have non-zero relevant terms corresponding 

to: 1^, #^ and 

The result are shown in the table 5: the second and third columns refer, respec- 
tively, to the fiducial values for the parameters and their 1 — a errors calculated from the 
Fisher matrix, while in the third column we report the shifts on parameters in the form of 
dOa/dfisoi that is the shift on the parameters per unit shift in fiso- The parameters more 
affected (labeled with ► in the table) are the scalar amplitude and spectral index Ag, Ug, 
the Hubble constant, parameterized by h, the cold dark matter physical density iOc and, as 
a consequence, the value of the sound horizon at radiation drag rs{Zfi). The fiducial value 
for the amount of isocurvature used in the Fisher calculation is fiso = —0.01 so, in the 
case that the universe were adiabatic, we would have a shift on parameters corresponding 
to Sfiso = +0.01. In this case it is interesting to compare the effect on the parameters 
6a induced by such a (small) shift to the expected error on the parameters: the values of 
this ratio (%) arc shown in the last column of table 5. The shift on the parameters, for a 
^fiso = 0.01, turns out to be bigger than the error (~ 120% of aa) for Vgizd), h and ojc and 
even more for Ug (— 150%). The amount of the shift in the cosmological parameters due 
to the presence of an isocurvature fraction found with this Fisher analysis is in agreement 
with what we found with the previous MCMC analysis, even if the base parameter set is 
different in the two applications (in the MCMC approach r and dng/dlnk were not left 
as parameters). In particular, we find in both case the same degeneracy direction for the 
parameters Hq and rs{zd) with fiso- 

While the shift is large, we expect the amplitude of this effect to be specific of 
the model we have considered {uiso = riad) and thus the amplitude of this shift to depend 
strongly on the choice adopted for the relation between the two spectral slope indices. The 
quantitative effect on As can also be affected by this choice. However, the effect on the 
other parameters such as rs{zd), is expected to be less sensitive, at least qualitatively, to 
this choice as the signal comes from a more localized region in multipoles I (see discussion in 
[34]). Thus while the effect on Us may be relevant to analyses geared towards determining 
the shape of the primordial power spectrum and investigating the implications for inflation, 
here we concentrate on the implications of a systematic shift on rs{zd) and parameters 
yielding the Universe's expansion history. 

4. Isocurvature and BAO: Future surveys forecast and implication for DE 

In light of the findings obtained so far, the aim of this section is to explore the effect 
on the BAO observables of allowing for an isocurvature fraction. The BAO observables 
are Eqs. (A. 6) and (A. 7): B\Az) = H{z) rs{zd) and B±{z) = D a{z) / r s{zd) ■ These are 
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Table 5: Fisher matrix for a Planck- like set-up - Shift on parameters due to the fact that /jgo is 
fixed to an incorrect value {6fiso 7^ 0). The fiducial value adopted here is fiso = —0.01. 





r IQ. VdiUc 


r ibiit^r 61 1 or ^ 


CJtiiff 

OlllH -jr-p 


| — l<5/i,o=0.oU/oj 


r 


U.Ui 


u.uzzyy 


U.iODo 


D. / y 


► Ug 


u.yoo 




U.DDiZ 


lOZ. / 


dus/dln k 


0.0 


0.0053 


-0.2292 


43.24 




0.840 


0.0077 


-0.0196 


0.25 


Ub 


0.02273 


0.000127 


-0.0005 


3.93 


► UJc 


0.1099 


0.00115 


-0.139 


120.55 


► h 


0.72 


0.00581 


0.71467 


123.0 


► As 


0.8169 


0.00825 


0.35752 


43.33 


fiso 


-0.01 


0.01123 


/ 


/ 


► rs(2;d) Mpc 


149.641 


0.359 


-44.215765 


123.19 



the quantities, respectively parallel and transverse to the line-of-sight, that future galaxy 
surveys like BOSS [59] and EUCLID [60] plan to measure with high accuracy. We quantify 
the effect that the presence of isocurvature could introduce on the expected values and 
compare it with forecasted errors on such quantities. In order to do so we use the code and 
the method described in [44]. The experimental specifications for both surveys are listed 
in table 3.2. We refer to section 3.1 for the calculation of the expansion history H{z) and 
the sound horizon at radiation drag rs{zci)- Similarly, we compute the co- moving angular 
diameter distance Da{z) from Eq. A.l. 

The results are summarized in Fig. 2 and Fig. 3. On the left panel we show the 
redshift dependence of the observable quantity {rs{zii)H{z) and Da{z)/rs{Z(i) respectively) 
and the forecasted error-bars for forthcoming (BOSS) and planned (Euclid) surveys. On 
the right panel wc show the % shift of the quantity (thin lines) and the forecasted errors 
(dashed thick lines, black and purple). We should distinguish two effects. 

The first is a systematic effect arising from the shift on the CMB-based rs{zd) deter- 
mination. This is shown as the red solid line at the 1% level for a A/iso = 0.08. This 
effect is negligible for a survey like BOSS but it is at the level of statistical errors for the 
Euclid H{z) over 7 or 8 redshift bins and at or above Euclid statistical errors on Da{z) 
over more than 10 bins. This indicates that the presence of an isocurvature component 
at a level well within the allowed range, if neglected, can introduce a systematic error in 
the interpretation of the BAO signal that is comparable if not larger than the statistical 
errors. Similarly, an isocurvature component, fiso, of a magnitude allowed at 3 — o"-level 
by Planck data could yield a 1 — cr systematic effect on B\\ and therefore H{z) and an fiso 
of a magnitude allowed at 2.5 — cr-level could yield a 1 — a systematic shift on B± and thus 
Da{z). 

One should keep in mind that a 1% shift on H{z) at z ~ 1 implies a shift of 0.05 on 
Wo and that a 1% shift on Da at 2; ~ 1 implies a shift of 0.07 on wq. 
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Figure 2: BAO observable H{z)rs{zd) and isocurvature (fiducial curvaton model with mixed 
ad+cdi modes fiso = —0.08 ): LEFT panel: the dashed orange line shows H{z) rs{zd) as a function 
of the redshift z for the fiducial curvaton model, the blue line for the fiducial curvaton model fitted 
with an adiabatic model (Fid. ISO-fit. AD) for a standard ACDM cosmology and the dot-dashed line 
for the same model but for a cosmology with varying dark energy equation of state w [the black 
line refers to the curve obtained by using the best fit parameters while the green curve by using 
the mean values for the parameters]. The purple error bars are the forecasted 1 — a errors for the 
BOSS survey [59], while the black ones are for the EUCLID survey [60]. The RIGHT panel plot 
shows the percentile variation of H{z)rs{zd) calculated by comparing the value of this observable 
as a function of z for the underlying fiducial curvaton model with the fiducial curvaton model fitted 
with an adiabatic model (Fid. ISO-fit. AD) for a standard ACDM cosmology (blue line) and with 
the same model but with varying w (dot-dashed lines). The dashed lines refer to the percent 1 — a 
errors on H{z) rs{zd) for the BOSS survey (purple) and EUCLID (black). The red straight line 
refers to the percentile variation of H{z) rs{zd) when only rs{zd) is shifted because of the presence 
of the isocurvature (Fid.ISO vs Fid.ISO-fit.AD ACDM, see section 3.1.1). 



The second effect is a mis-match of the expansion history as inferred from CMB and 
as measured by BAO surveys. The thin lines in Figs. 2, 3 show the CMB-predicted BAO 
observables (which combine the shift on the expansion history and the shift on the sound 
horizon at radiation drag). The solid blue line corresponds to the Fid.ISO-fit.AD case. 
The dot-dashed lines are the best fit (black) and mean value (green) for the Fid. ISO- 
fit. ADw case. The mis-match is clearly above the BAO surveys statistical errors over 
a wide range in redshift for the transversal measurement for BOSS and for both, radial 
and transversal measurements for Euclid. The mis-match between CMB predicted and 
the measured expansion histories has been proposed as a signature for deviations from a 
ACDM cosmology in the form of deviations from Einstein's gravity e.g., [62, 63], couplings 
in the dark sector e.g., [61] or time evolving dark energy. Here we add another possibility 
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Figure 3: BAO observable l and isocurvature (mixed ad+cdi model): LEFT panel: the 

dashed orange line shows i^Alf) as a function of the redshift z for the fiducial curvaton model, 
the blue line for the fiducial curvaton model fitted with an adiabatic model (Fid. ISO-fit. AD) for a 
standard ACDM cosmology and the dot-dashed line for the same model but for a cosmology with 
varying dark energy equation of state w [the black line refers to the curve obtained by using the 
best fit parameters while the green curve by using the mean values for the parameters] . The purple 
error bars are the forecasted 1 — a errors for the BOSS survey [59], while the black ones are for the 
EUCLID survey [60]. The RIGHT panel plot shows the percentile variation of calculated by 

comparing the value of this observable as a function of z for the underlying fiducial curvaton model 
with the fiducial curvaton model fitted with an adiabatic model (Fid. ISO-fit. AD) for a standard 
ACDM cosmology (blue line) and with the same model but with varying w (dot-dashed lines) where 
the black line refers to the curve obtained by using the best fit parameters while the green curve by 
using the mean values for the parameters. The dashed lines refer to the percent 1 — a errors on 
for the BOSS survey (purple) and EUCLID (black). The red straight line refers to the percentile 
variation of DA{z)/rs{zd) when only rs^zj) is shifted because of the presence of the isocurvature 
(Fid.ISO vs Fid.ISO-fit.AD ACDM, see section 3.1.1). 

to this list: an incorrect assumption about the nature of primordial perturbations. 

We should note however that including fiso as a parameter in the CMB analysis solves 
this issue. Should there be a small isocurvature contribution, adding the fiso parameter to 
the analysis vi^ill recover correctly and the expansion history vi^hile the error-bars degrade 
by quarter of a sigma. 

5. Conclusions 

In this paper we have investigated how a biased determination of of the sound horizon 
due to a incorrect assumption about the nature of the initial perturbations affects cosmo- 
logical parameters measurements from future BAO data. It is important to keep in mind 
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that analyzing the CMB data with the prior assumption of purely adiabatic initial con- 
ditions when the real universe contains even a small isocurvature contribution, could lead 
to an incorrect determination of the cosmological parameters. In fact, the presence of an 
isocurvature component changes the shape and the location of the CMB acoustic peaks, 
mimicking the effect of parameters such as ^mh^ , Hq and w. Here concentrated on wether 
such an effect has an impact on standard rulers, like the sound horizon at recombination, 
inferred from CMB observations. This issue turns out to be very delicate and important 
especially in view of the next generation of galaxy surveys which aims at probing with high 
accuracy the late time expansion and thus the nature of dark energy by means of BAO at 
low redshift {z <2). 

We find that the presence of an isocurvature contribution, of a magnitude still allowed 
by present data, can affect both the size of the sound horizon at radiation drag and the 
CMB-inferred expansion history, in a non- negligible way. For our forecasts we used an 
experiment with the characteristics of the Planck mission. We find that the systematic 
error introduced in the sound horizon at radiation drag propagates into a systematic error 
on quantities like H{z) and Da{z) that can be comparable or larger than the statistical 
error over a wide range of redshift bins for future BAO experiments. 

A crucial assumption of the BAO technique is the possibility to measure accurately 
and robustly the acoustic sound horizon at radiation drag. Generally CMB observations 
yield a robust determination of such standard ruler although there could be possible sys- 
tematic effects introduced by deviations from the standard evolution of the early universe. 
In this paper we have presented another possible deviation from the minimal, standard cos- 
mological model that could affect the BAO interpretation as standard rulers if not taken 
into account: the presence of an isocurvature contribution to the initial conditions for 
primordial perturbations. 

In addition the neglected presence of an isocurvature component in the initial condi- 
tions introduces a mis-match between the expansion history as inferred from CMB and 
as measured by BAO or Supernovae surveys. The mis-match is above the BAO surveys 
statistical errors over a wide range in redshift for both on-going and future surveys. Such a 
high redshift-low redshift mismatch between CMB predicted and the measured expansion 
histories has been proposed as a signature for deviations from a ACDM cosmology in the 
form of deviations from Einsteins gravity, couplings in the dark sector, or time evolving 
dark energy. Here we add another possibility to this list, an incorrect assumption about 
the nature of primordial perturbations. 

In the paper we have used the curvaton scenario as a working example for a model that 
gives rise to a small fraction of correlated isocurvature perturbations. The curvaton has 
Uad = niso-, therefore our findings are quantitative only for this case, but hold qualitatively 
for other models. 

In order to recover an unbiased determination for the sound horizon and dark energy 
parameters, a component of isocurvature perturbations must be included in the model when 
analyzing CMB data. Fortunately, doing so does not increase parameter errors significantly. 
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A. BAO and Iso curvature 

The baryon acoustic oscillations (BAO) responsible for the acoustic peaks structure im- 
printed in the CMB are predicted to be present in the late-time clustering of galaxies as a 
series of weak modulations in the amplitude of fluctuations as a function of scale [44] . 

In general, the line-of-sight (ry) and transverse (rj_) co-moving sizes of an object or a 
feature at redshift z are related to the correspondent observed sizes Az (redshift slice) and 
A0 (angular separation) by the expansion history il[z) and the angular diameter distance 
through: ry = and r_L = D^(2;)A0, being the co-moving angular diameter 
distance: 

r dz' 

^.(.)^.^^ (A.1) 

and 

H{z) = Ho ^^0,(1 + 2)3(1+-^) (A.2) 

where Hq is the Hubble parameter at the present time and fij and Wi are, respectively, the 
density parameter and the equation of state parameter of a given species i. 

The key parameter associated to the BAO in the primordial photon-baryon fluid, which 
leave an imprint in the large scale matter distribution, is the size of the sound horizon at the 
radiation drag redshift Zd, i.e. at the time when baryons were realized from the photons. 

with the sound speed: 

c,{z) = '^^^ ^ , (A.4) 

being Oj, and Q^, respectively, the baryon and photon density parameters at the present 
time. The function E{z) = (^^ ni(l-|-z)3(i+"'»)) 2 up to the time of radiation drag accounts 
for the contribution of radiation (photons and neutrinos) and matter (baryons and cold 
dark matter) so that takes the form: 



E{z) ~ V^r{l + z)^ + n^{l + z)\ (A.5) 

where $7^ = ^^7Q(l + 0.2271iVcfj) is the radiation term, $7^,, is the present photon density 
parameter and A'^eff the neutrino effective number that we set to the standard choice of 3 
neutrino species. 
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The CMB could in principle provide such a 'standard ruler', so that, by measuring the 
angle subtended by this ruler as a function of redshift we can map out the angular diameter 
distance Da{z) and by measuring the redshift interval associated with this distance we can 
map out the Hubble parameter H{z). This is the key idea on which the next generation 
of future BAO galaxy surveys is based. 

We will thus be interested in the observables, respectively along and perpendicular to 
the observer line-of-sight: 

B^^{z) = H{z)rsizd) (A.6) 

and 

In view of high precision data it is mandatory to test how robust is the method and the 

determination of the standard ruler. In particular it will be important to know/account 
for all the possible systematics that could affect such measurement. As already pointed 
out in the Introduction, relaxing the hypothesis of adiabaticity for the initial conditions 
compromises the precision of parameter estimation: for example, in the case of the sound 
horizon at radiation drag Vs^za), the presence of an isocurvature contribution does not 
change the dependence on the cosmological parameters, but it will produce a difference in 
the measured value with respect to the one obtained assuming a purely adiabatic model. 
Analogously the presence of an isocurvature contribution would produce an "adjustment" 
in the cosmological parameters which could lead to a wrong interpretation of the expansion 
history parameter. 

B. Fisher matrix for CMB 

The Fisher matrix is defined as [54]: 



dOidOj 

where C is the (gaussian) likelihood for a set of parameters 0: 



~ (27r)'^/2|detCov|V2''''P 



2 



(B.2) 



D and y = y{9) are respectively the data set and the fiducial theoretical model and Gov 
is the covariance matrix: CoVjj = {{Di — yi){Dj — yj)). 

In general the estimate of the covariance for the parameters is given by: cr?^- > {F~^)ij 
and the marginalized error (i.e. calculated with the full covariance matrix) on a given 
parameter 6i is: 

ae, > ViF-'h- (B.3) 

So, once computed the covariance matrix and having a fiducial model with a known de- 
pendence on the parameters, the Fisher matrix gives the expected errors. This is a very 
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powerful tool extremely useful in experiments design which enables to forecast the param- 
eter errors before doing the experiment. 

In the case of a full sky and noiseless CMB experiment the gaussianly distributed 
signal in the sky can be described by the vector {aj^, af^, af^), with a^, X = T,E,B the 
spherical harmonic coefficients for the temperature, the E and the B polarization modes. 
In this case the Fisher matrix will take the form: 

XY I ^ ^ 

where Q is the CMB angular power spectrum: {aemO'£im') — ^U'Smm'Ce and Aif^, X,Y = 
TT, EE, TE, BB, ... are the elements of the matrix: 



Me 



21 + 1 



( {Cj^f (C™)2 cj^cf^ \ 

Cj^CP Cf'Cj^ l/2[(Cf + CrcJ^] 

V {cf^fj 



which takes into account the correlation between TT, EE and TE. 

C. Isocurvature modes and the CMB: parameter degeneratios 

The main effect of pure isocurvature perturbations on the CMB is a | phase shift on the 
acoustic oscillation with respect to the adiabatic case which moves the peak structure in the 
CMB temperature spectrum to smaller scales. This could be qualitatively described as a 
phenomena of 'compensation' [48]. Any source of isocurvature density perturbations must 
be compensated to keep the total density fluctuations small, this means that there must 
an anticorrelation between the isocurvature source and the photon density perturbations 
(being the photons the dominant dynamical component at this stage). 

More specifically, in the adiabatic case the photon-baryon fluid begins to compress itself 
due to its self-gravity. As pressure tries to stop the compression, the potential decays. 
The fluid is left in a highly compressed state and self-gravity acts to enhance the first 
compression of a cos(ks) series. Peaks occur at: 

ei'' = meA m = 0,1,2,... (C.l) 

where £a is the ratio between the co-moving angular diameter distance to the last scattering 
surface Da{z*) and the sound horizon at last scattering rs{z*) (Eqs. (A.l) and (A.3)): 

U-.^. (C.2) 

rs{z*) 

Baryon drag enhances all odd peaks. 

In the isocurvature case the potential fluctuations vanish initially and then grow in 
anticorrelation with the photon fluctuations until horizon crossing. Because of photon 
pressure near horizon crossing the compensation mechanism that keeps the curvature van- 
ishing breaks letting the potential fluctuation grow from zero. To compensate this the 



-19- 



photons rarefy inside the potential wells. At horizon crossing the photon pressure resists 
the accompanying rarefaction and the fluid turns around and begins falling into the po- 
tential wells. The self-gravity of the photon-baryon fluid drives the sine rather than the 
cosine oscillation. The position of the acoustic peaks is given by: 

^^ = {^-\)^A (C.3) 

and baryon drag enhances all even peaks. 

Therefore, accounting for an isocurvature contribution induces a shift in the position 
of the CMB peaks to smaller scales: the more isocurvature fraction, bigger the peak dis- 
placement. However, this effect can be compensated by a change in other cosmological 
parameters, namely the Hubble parameter Hq. Qualitatively an enhancement in yields 
in fact a shift of the peaks position towards larger scales while a decrement results to a 
shift to smaller scales which can thus mimic or compensate the presence of isocurvature 
modes. For this reason the two parameters jiso and B.q are said to be degenerate. For 
a more detailed discussion on isocurvature and degenerations with the CMB parameters, 
in particular for the more complicated cases of the mixture of adiabatic and isocurvature 
modes, see for example: [37, 34, 41]. 
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